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PREFACE 


The purpose of this report is to facilitate fracture control in glass and ceramic systems. This report 
is intended for those materials engineers who are not thoroughly familiar with fatigue testing in glasses 
and ceramics but who need to develop design-allowable estimates for structures of these materials. 
Since this type of testing is relatively new and different from the testing in metal systems, a rationale 
for the test methodology is presented contrasting it with metal systems, which are more familiar to 
many materials engineers. The procedures developed are applicable to materials with natural process- 
ing and fabrication flaws. Although the procedures can be applied to materials with controlled flaws, 
no attempt is made at addressing advances made in the area of controlled flaw (indentation) testing. 
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Introduction 

Determining the safe lifetime of a glass or ceramic structure using fracture mechanics has been 
the subject of intensive investigation since its successful use in ensuring the reliability of aircraft and 
spacecraft windows (reference 1). Theory states that the independent strength-related property of these 
materials is their critical stress-intensity factor and that this material constant is proportional to the 
product of the applied stress at failure, times the square root of the failure-initiating crack size. Thus, 
these brittle materials do not possess a single intrinsic strength, but instead, they possess a strength 
distribution that depends on the distribution of the largest flaws in the material. Since these materials 
typically experience subcritical crack growth when stressed in tension in the presence of moisture, 
predicting a safe lifetime for a piece of hardware requires a determination of crack dynamics in the 
use environment as well as a determination of the largest preexisting crack size in critical areas of 
the hardware. Therefore, predicting a safe lifetime for these materials requires an analysis of the follow- 
ing three factors: (1) the size of the largest crack in the areas of highest stress in a piece of hardware, 
(2) how fast the crack grows under stress, and (3) how big the crack can get before it causes failure. 

Unfortunately, the first two factors must usually be arrived at indirectly. In contrast to metals, 
determination of the largest crack in a piece of glass or ceramic is complicated because by the time 
a crack is observable in these materials, using most presently available nondestructive test techniques, 
the material is useless as a load-bearing member (Appendix A). Thus, the approach usually followed 
is to proof-test the hardware, thereby indirectly ensuring that cracks are less than a certain size. This 
approach assumes that the expression for the material’s critical stress-intensity factor is valid and that 
the cracks do not grow during proof-test unloading. The latter is usually ensured by running the test 
in a dry environment and unloading rapidly (reference 2). 

The next step in lifetime prediction is to determine crack dynamics. Large-crack propagation tech- 
niques, so useful for metals, are typically done first. These tests are quick, conservative of material, 
and not labor intensive. However, their use is complicated by some evidence that the crack propaga- 
tion constants determined in this manner are not always the same as those derived indirectly by measur- 
ing the strength of specimens with small cracks (reference 3). Therefore, crack velocity measuring 
techniques such as double cantilever beam and double torsion (reference 4) are not always applicable, 
especially to materials with a polycrystalline microstructure. To circumvent this difficulty, indirect 
statistical techniques have been developed. These techniques are stress-rate testing (dynamic fatigue) 
and time-to-failure testing (static fatigue). These tests are usually modulus of rupture measurements, 
although biaxial tension measurements have some advantages (reference 5). 

In a hardware design effort, time is usually a limiting factor. The usual progression is to begin 
with large-crack propagation techniques and proceed to dynamic fatigue and finally to static fatigue 
as time permits. The end result of these procedures is to produce a design diagram. 



A design diagram is a parametric plot of time to failure versus applied stress. In using these three 
techniques to produce design diagrams, time-to-failure predictions are extremely sensitive to experimen- 
tal uncertainty in the fatigue parameters. Statistical analyses have been developed in the literature 
(references 3, 6, 7, and 8) to assess this uncertainty so that confidence intervals can be applied to 
time-to-failure estimates in design diagrams. Since this report is written from the point of view of 
a materials engineer concerned with developing design-allowable estimates, a good deal of emphasis 
is placed on confidence interval estimates. 

This report illustrates how to implement fatigue data to graph design diagrams with error bars 
using an APL computer program. The error analyses used are called the median analysis and the 
homologous stress-ratio analysis for dynamic and static fatigue (reference 3). The report uses, as an 
example, testing that was performed on a machinable glass-ceramic* used structurally in the Energetic 
Gamma-Ray Experiment Telescope (EGRET) on the Gamma-Ray Observatory (GRO), which is 
scheduled to be launched in 1988 (references 9 and 10). 

Discussion 

A design diagram is a parametric plot of lifetime versus applied stress. To arrive at a design diagram, a starting 
point is the expression for a stress-intensity factor in opening mode I, which is the tensile mode,^ 

Kj = Yaa Vl (1) 

where Y is a geometric constant, taken as \fi f for surface cracks (reference 11), and a is the far-field stress'^ applied 
to a crack of linear dimension, a. Mode I stress intensity is typically the only one considered because glasses and ceramics 
do not undergo much plastic deformation, and consequently they fail in tension rather than in shear. 

In engineering terms, the theory states that these materials break when the product of a and a 1/2 exceeds a critical 
value, K JC (reference 12). Thus, the independent strength-related property of these materials is their critical stress in- 
tensity factor given by the equation 


K. r = Ya a 1/2 

IC c c 

where a and a are the stress and crack size at failure. 

C C 

According to this theory, a material does not possess a single characteristic strength. It possesses a distribution 
of strengths that reflects its distribution of cracks or flaws (Appendix A). If the material does not experience sub- 
critical crack growth, measurement of K JC is where materials strength testing ends and therefore a design diagram 
is not needed. The part is designed, processed, and maintained in service so that the biggest crack in it is small enough 
to give the required strength. If, however, the material does experience subcritical crack growth under stress, then 
an estimate must be made as to how long it takes the largest crack in a piece of the material to grow to critical size. 
In estimating the time required for a crack to grow subcritically, a simple power law relation between subcritical crack 
velocity, V, and stress intensity, K p is most often assumed (reference 2). The relation is 

V = AK t N (3) 

where A and N are material/environmental constants. When this is done, it has been shown (Appendix B) that frac- 
ture stress, S, at a constant stressing rate, a, is given by the expression 

S N+l = B(N+ l)Sf- 2 a (4) 


♦Macor™, Corning 9658 glass-ceramic 

tin contrast to sliding shear, mode II and tearing shear, mode III (reference 1 1) 

tAs distinguished from the concentrated stress at the root of the crack tip (references 12 and 13) 
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where 


B = 2/ [A Y 7 (N-2)K£ 2 ] 
S. = inert strength* 


Taking the natural logarithum of equation (4) and dividing by N+ 1 linearizes the equation giving the expression 

In S = [In B +ln {N+ 1) + (N- 2) In S. + In 6] (5) 

For a given inert strength, the increase in observed strength, S, as a function of stressing rate, predicted by equation 
(3) is called dynamic fatigue. ^ 

Similarly, it has been shown (Appendix C) that time to failure, t f , under a constant applied stress, a g , is given 
by the expression 

t. = BS N - 2 a N (6) 

/ ' ° 

Linearizing equation (6) by taking its natural logarithum gives the expression 

In / = In B + (N-2) In S. - Nln cr a (7) 

For a given inert strength, the delayed failure behavior predicted by equation (6) is called static fatigue, and 
phenomenologically it is the same as stress rupture in metals. 

The fatigue parameters N and B can be directly determined by measuring K [C (which is needed to calculate B) and 
by using equation (3) to analyze large-crack velocity data (reference 14). Alternatively, since 


S. 

I 


= K ic /Y °i 


( 8 ) 


where a. is the initial size of the failure-initiating flaw, these constants could be measured directly from fatigue data 
if there were some way of measuring, on the same specimen, both the initial size of the preexisting flaw that causes 
failure and the fatigue behavior (i.e., strength for a given stressing rate (equation 5) or time to failure for a given 
applied stress (equation 7)). Unfortunately, nondestructive test techniques in ceramics have not been developed to 
the point that the initial, preexisting, flaw size that causes failure can be measured consistently. One way of measuring 
that flaw size, however, is to break the specimen in such a manner that no subcritical crack growth occurs (i.e., measure 
its inert strength, S). The failure-initiating flaw size is then determined using equation (8). 

This technique has led to a determination of the fatigue parameters A and B using equations (5) and (7) by substituting 
distributions for the variables and assuming that the distributions are equivalent and partitioned in the same way 
(references 15 and 16). Thus, for example, it is assumed that median values (S. , t f or S j , §) of a test sample of 
nominally identical specimens can be used for S ( - , tj , and S in equations (5) and (7) to arrive at values for N and 
B. The median value is used since it is particularly suited to tests in which only a fraction of the total number of 
specimens is broken, such as static fatigue where it is often impractical to wait for all the specimens to break. 

Thus, to determine the fatigue parameters N and B using fatigue data, a population of nominally identical test 
specimens is made. The specimen population is then divided into test samples containing equal numbers of specimens. 
For a full complement of fatigue tests, this laboratory uses a population of 210 specimens divided into 7 sample sets 
of 30 specimens each. First, a set of inert strengths are measured, and then the remaining samples are used to measure 
the median strengths at three stressing rates for dynamic fatigue and three median times to failure at three applied- 
stress levels for static fatigue. 


♦Inert strength is the term given to the strength of a ceramic or glass material when there is no subcritical crack growth. This strength is usually 
measured so quickly that subcritical crack growth does not have a chance to proceed. Otherwise, it is measured in an environment that pre- 
cludes crack growth (i.e., a very dry or a very cold environment). 

^This is in contrast to metals in which the term fatigue is usually reserved for a cyclical process in which a brittle, cold worked zone is created 
at the crack tip and carried in front of the crack until a critical size is reached at which point fast fracture occurs. 
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For dynamic fatigue, a least squares plot of In S versus In a is made using the median inert strength S. in equation 
(5). Similarly, for static fatigue, median times to failure are measured at various stress levels, and a least squares plot 
is made of equation (7) using median values. From these regression analyses, the fatigue parameters N and B can 
be evaluated from the slopes and intercepts, respectively. With acquisition of the fatigue parameters N and B, equa- 
tion (7) can be used to create a design diagram. Using this equation, a plot of In t f versus In a g is a series of straight 
lines with a slope -N, which are parametric in S. or equivalently in initial crack size a . , using 0 equation (8). Neither 
of these plots is very useful since the parameters S. and a. are not readily determined in a piece of hardware to verify 
the initial quality of the piece. 

However, the inert-strength distribution of glasses and ceramics can usually be approximated by a two-parameter 
Weibull distribution 



where S g and m. are called the scaling parameter and the Weibull modulus, respectively, and Fis the cumulative prob- 
ability of failure. To evaluate the Weibull parameters m. and S g , a least squares plot is made of the logarithmic restate- 
ment of equation (9): 

x F = m. In S. - m. In S o (10) 

where x p = In In (1/1-F). Equation (7) can then be restated as a design diagram, parametric in probability of failure, 
by solving equation (10) for In S j and substituting it into equation (7), which gives the expression 

In t f = In B + [x F + m. In SJ - N In a g (11) 

This is the form of the design diagram that has been found most useful for initiating design. Since the equation does 
not take into account differences in the stress distribution between test specimens and the hardware, or size differences, 
it is truly an initial step in an evolving and iterative process as discussed in Appendix D. 

Although the median value technique for analyzing fatigue data is quite straightforward to apply, it makes ineffi- 
cient use of the data since the least squares analysis uses only one value per sample of 20 to 30 specimens. Recognizing 
the need for more efficient utilization of fatigue data, several researchers (references 17 and 18) have suggested a method 
of data reduction that is based on homologous stress ratios a HS and a HD (reference 19) defined by 


a HS £ an d a HD £ 

i i 

(12) 

Rewriting equations (5) and (7) in terms of these ratios gives 


ln °hd ~ a/+i [ ln B + 1" (N+l) + In ] 

(13) 

and 


In (t^. Sf) = InB- Nlna HS 

(14) 


By ranking the S data from weakest to strongest for a given a and then by pairing these data with equally ranked 
S. data, a least squares analysis of In o HD versus In (<j/S?) will give N and B from the slope and intercept, respective- 
ly. Similarly, from a least squares analysis of In (t f S 2 ) versus In a HS , the parameters N and B are determined from 
the slope and intercept, respectively, of the static fatigue data. Note that times to failure are ranked from shortest 
to longest and are paired with increasingly ranked inert strengths from a sample of specimens whose number equals 
the number of specimens under load in the static fatigue test. If it is necessary, starting sample sizes are made equal 
by using a random number generator to discard the excess number of specimens. 

Since all the data are used rather than just one value per sample of 20 to 30 specimens, this method increases the 
confidence in N and B as compared with the median value technique (reference 7). The main advantage of using a 
homologous ratio analysis, however, has been that it points out anomalous behavior, particularly for static fatigue 
data. An example of this appears in the section discussing the use of the program. 
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Time-to-failure predictions using equation (11) are very sensitive to uncertainty in the size of the failure-initiating 
crack and uncertainty in how fast that crack grows. This is reflected in experimental uncertainty in S ( . and the fatigue 
constants TV and B. Expressions for this experimental uncertainty and its effect on the uncertainty of In iy can be resolved 
(references 3 and 6). Determining these expressions entails application of the error-propagation rule (references 3, 
6, 20, and 21) to equation (7). The error-propagation rule gives the uncertainty of a quantity in terms of the uncertain- 
ty of the variables used to calculate that quantity. In engineering terms, the rule states that if x is a quantity that 
is a function of at least two other variables n and v (i.e., x = fin, v , ...)), then the uncertainty of x in terms of its 
variance V(x) is given in terms of the variance and covariance of the variables n and v by the expression 

F(x) = FOt) 0 2 + nv) 0 2 + 2 Cov + - • ■ ’ (15) 


there being additional terms and cross terms for additional variables. Since the standard deviation, a, of x is the square 
root of the variance of x, a plot of x with 95 percent reproducibility limits would be the value of x ± 1 .96 times the 
standard deviation of x. Similarly, the 90 percent reproducibility limits would be ± 1.645 standard deviations, and 
so on (reference 3). 

Application of the error-propagation rule to In /y in equation (7) gives the expression 



(16) 


Since a is supposed to be known exactly, terms involving the variance and covariance of In o a have been eliminated 
from equation (16). Furthermore, since TV and B are functionally unrelated to S., covariance terms between S. and 
those variables are also zero. Evaluating the partial derivatives in equation (16) results in the expression 


F(ln t f ) = F(ln B) + (ln-L) V(N) + (TV-2) 2 F(ln S.) 




+ 2 Cov (TV, In B) In & 


(17) 


In fatigue testing, expressions for F (In B), V(N), Cov (TV, In B), and F (In S.) can be found in a number of ways, 
depending on the assumptions made concerning the type of distribution describing jy, S, and S., and on the type of 
analysis used (i.e., median or homologous stress-ratio analysis). An expression for F(ln S ( .) is derived in Appendix 
E, assuming that S. is described by a two-parameter Weibull distribution. Alternatively, F(ln S ( .) can be calculated 
directly from the measured inert-strength data using the commonly accepted expression for variance (reference 20) 

F(ln S.) = a 2 = — - £ [(In S.)j - (hT^)] 2 (18) 

T=i 

where In S ( . is the mean value for In S. (i.e., l/n (In S.)p. Either expression for F (In S ( .) can be used in the median 
strength analysis or the homologous stress-ratio analysis. The other three terms in equation (17) can be calculated 
in two ways for the median analysis, but only the first of these ways can be used to calculate these terms for the 
homologous stress-ratio analysis as shown below. 

Since both the median and homologous ratio techniques deal with linear equations, the method for arriving at 
expressions for F(ln B), F(TV), and Cov (TV, In B) can be generalized in terms of the usual equation for a straight line: 

y = ax + b (19) 

where y and x are the dependent and independent variables, respectively, while a and b are the slope and intercept. 
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In linear regression or least squares analysis, it is possible to calculate the variance of the least squares slope and 
intercept parameters as well as their covariance by using the /-paired values of the dependent and independent variables 
using the expressions (references 3, 20, 21, and Appendix F): 


F(C) JR (x) 

(20) 

V(b) = V(a) [R(x) + x 2 ] 

(21) 

Cov (a, b) = -V(a) x 

(22) 


where j 

° 2 = JJi ?, + 

R(x) =-j£ {x.-xf 

i=i 

_ 1 ' 

X = — L x ' 

/= i 

In the first method for arriving at expressions for F(ln B), V(N), and Cov(N, In B), the error-propagation rule 
is applied to expressions for a and b after they have been solved for N and In B, respectively. This gives equations 
for V(N) and F(ln B ) in terms of F(c) and V(b). The covariance of N and In B is then evaluated in terms of Cov(a, 
b) by using the differential expression for covariance found in reference 3 and developed in Appendix F. This process 
is performed in Appendix G for median analysis and in Appendix H for the homologous stress-ratio analysis. 

The second method for arriving at expressions for F(ln B), V(N), and Cov(N, In B ) applies only to the median 
analysis. The method assumes that t f , S, and S. are all described by two-parameter Weibull distributions. Making 
that assumption, it is possible to derive expressions for the parent sample variance, a 2 , for the median fracture strength 
S, and the median time to failure^., in terms of the number of specimens tested, the Weibull modulus, and the fatigue 
parameter, N (reference 3). In Appendix I, these expressions for F(ln B), V(N), and Cov(N, In B) are reproduced 
from reference 3 for convenience. 

The first method of error-bar estimation (i.e., using only measured data) is preferred for median analysis rather 
than the second method (i.e., assuming Weibull distributions, for reasons outlined in Appendix E). 

Annotated Programs 

The analysis just reviewed has been incorporated into an APL computer program called DESIGN. The program 
uses two functions available in APL public libraries and three functions developed to use with DESIGN as subroutines. 
These functions are L TF, SVLTF, and DVLTF. The public library functions are PCF and GRAPH, which are found 
in reference 22. 

The function PCF is a least squares routine that is used for determining the slope and intercept of the linear equa- 
tions used to evaluate Weibull parameters and fatigue constants. In the program DESIGN, the output of PCF is a 
two-element vector: The first element is a slope, and the second an intercept. 

The function GRAPH is used to plot the design diagrams. Two types of diagrams are developed by DESIGN: One 
diagram is a plot of In t f versus In o a for various probabilities of failure without error bars. The other is a plot of 
In i versus In a g for a single probability of failure with error bars. 

In the program DESIGN, nemonic abbreviations are used wherever possible. The letter “F” after a designated 
variable such as applied stress “AS,” time to failure “TF,” log-applied stress “LAS,” failure stress “FS,” strength 
“S,” etc., indicates a vector. The letter “F” appearing before a variable such as log time to failure “LTF,” fatigue 
parameters “N” and “B,” or log inert strength “LSI” indicates the variance of that variable. The letters “S” and 
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“D” at the beginning of a function or variable indicates static or dynamic fatigue. Thus, VLTF means “variance 
log time to failure,” whereas SVLTF means “variance log time to failure using static fatigue data,” and DVLTF means 
the same thing using dynamic fatigue data. 

The commands SADESIGN~- 16, 93, 203; SADVLTF+-9-, and SASVLTF*-9 are stored in the work space named 
FATIGUE in order to regulate the program and allow data input. DESIGN is intended to be an interactive program 
in which the user can compare the results obtained by using different test methods and data reduction schemes. When 
the program stops before line 93, the results of the Weibull analysis are obtained by typing either WEIBULL or MI, 
SO. Similarly, when the program stops before line 203, the values of N and B are obtained by typing these letters 
or the words STATIC or DYNAMIC, whichever is appropriate. Values for V(N), F(ln B), F(ln S), and Cov(N, In B) 
are printed automatically during the execution of the functions SVLTF and DVLTF. 

A sequence of operations in the execution of DESIGN, SVLTF, DVLTF, and LTF is listed in Tables 1 through 
4. Following the tables, the programs are annotated in terms of the equations used in the text of this report. 


Table 1 

Operational Flow of DESIGN 


Step 

Operation Performed 

1 

Receive instruction about entering inert-strength data (lines 1-15). 

2 

Stop before line 16 to enter inert-strength data. 

3 

Calculate Weibull parameters m. and S g (lines 16-23). 

4 

Calculate variance log inert strength, not assuming a particular distribution, using measured 
inert-strength data (lines 24-25). 

5 

Receive instruction about entering fatigue data (lines 26-86). 

6 

Calculate median inert strength assuming a Weibull distribution (line 87). 

7 

Initialize x and y values for a least squares determination of slope and intercept (lines 88-92). 

8 

Stop before line 93 to enter fatigue data. 

9 

Calculate N and B using static fatigue median value data (lines 93-109). 

10 

Calculate N and B using static fatigue homologous stress-ratio data (lines 110-130). 

11 

Calculate N and B using dynamic fatigue median value data (lines 131-147). 

12 

Calculate N and B using dynamic fatigue homologous stress-ratio data (lines 148-166). 

13 

Receive instruction on formatting GRAPH to produce design diagrams (lines 167-196). 

14 

Receive instruction on producing a table of graphed data (lines 197-201). 

15 

Stop before line 203 to format GRAPH. 

16 

Graph design diagrams. 
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Table 2 

Operational Flow of SVLTF 


Step 

Operation Performed 

1 

Receive instruction (lines 1-8). 

2 

Stop before line 9 to decide the type of error analysis to be performed. 

3 

Calculate V(N), F(ln B ), F(ln S), and Cov(N, In B) assuming Weibull distributions and me- 
dian value analysis (lines 9-27). 

4 

Calculate V(N), F(ln B), and Cov(N, In B) using median value analysis but not assuming 
a particular distribution (lines 29-31). 

5 

Calculate V(N), F(ln B ), and Cov(N, In B ) using homologous ratio analysis (lines 33-35). 

6 

Print V(N), F(ln E), Cov(N, In B), and F(ln S). 

7 

Calculate variance log time to failure (line 46). 


Table 3 

Operational Flow of DVLTF 


Step 

Operation Performed 

1 

Receive instruction (lines 1-8). 

2 

Stop before line 9 to decide the type of error analysis to be performed. 

3 

Calculate V(N), F(ln B), F(ln S.), and Cov(N, In B) assuming Weibull distributions and me- 
dian value analysis (lines 9-29). 

4 

Calculate V(N), F(ln B), and Cov(N, In B) using median value analysis but not assuming 
a particular distribution (lines 31-34). 

5 

Calculate V(N), F(ln B), and Cov(N, In B) using homologous ratio analysis (lines 36-39). 

6 

Print V(N), F(ln B ), Cov(N, In B), and F(ln S). 

7 

Calculate variance log time to failure (line 50). 


Table 4 

Operational Flow of LTF 


Step 

Operation Performed 

1 

Calculate inert strength (line 1). 

2 

Calculate log time to failure (line 2). 
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VDESIGN[n]V 

VDESIGN;FV;Xl V; Yl V;J3;Xl ; WEIBULL;SIM; Y2 V;X2 V;X2;STA TIC; Y3 V;X3 V;X3;DYNAMIC;C90;LL;UL;X 
‘SI UNITS ARE USED FOR DATA ENTRY. STRENGTHS ARE IN’ 

‘PASCAL UNITS AND TIME IS IN SECONDS ’ 

‘ENTER INERT STRENGTH DATA STARTING ’ 

‘WITH THE LOWEST STRENGTH S( 1) AND ENDING WITH THE ’ 

‘ HIGHEST STRENGTH S(K) WHERE K IS THE NUMBER OF ’ 

‘SPECIMENS, FOR EXAMPLE ’ 

‘SI — 10.183F7’ 


‘S25 — 11 .596F7 ’ 

‘END THE DATA INPUT BY ASSIGNING A VALUE TO THE’ 

‘ VECTOR SV AS FOR EXAMPLE ’ 

‘SK— S1,S2,S3,...S25 ’ 

Value assigned to inert-strength vector SV = SI, S2,...S25 
‘FINALLY TYPE THE COMMAND —Cl TO CONTINUE THE PROGRAM' 


C\:FV— ((ipSV)-0.5) + pSV 

Value assigned to probability vector FV = (1—0.5 )/n, (2-0.5)/n, .... (n— 0.5)/n 
XIV— « SV 
XIV = In SV 
YIV— S ® 1 + 1 —FV 
YIV = In In (1/1 -FV) 

73 — 2, PSV 

J3 = a two-element vector : first element 2, second element PSV (25 in this example) 
XI— 73 PX1V.Y1V 

XI = a two-row, 25-column matrix {in this example ): first row XI V, second row yi V 
WEIBULL — 1 PCF XI 


YIV XIV 

Xp = m. In S. - m. In S Q , text equation (10) 
slope intercept 


WEIBULL = a two-element vector containing slope and intercept of text equation (10) 
MI— WEIBULU 1] 

MI = m f = first element of the WEIBULL vector 
SO — * — ( WEIBULL[2]) +MI 

£ = ^ -( 1/mp X ( intercept of text equation 10) 

SUM*— +/X1K 


SUM = E On S ). 

7-1 1 

VLSI— (+/{X\ V- SUM+ ( P S V)).2) + (( P S V) - 1) 

1 n 

VLSI = V(ln S.)=— x L [On S.).- (In S.)] 2 

•ENTER FATIGUE DATA: IF MEDIAN VALUE ANALYSIS IS BEING USED,’ 
•THEN THE DATA IS ENTERED IN ORDERED PAIRS BY EITHER ASSIGNING ’ 
' VALUES TO THE MEDIAN TIME TO FAILURE VECTOR, TFVM, AND ’ 

‘TO THE APPLIED STRESS VECTOR, ASV, FOR STATIC FATIGUE’ 

‘OR ELSE ASSIGNING VALUES TO THE MEDIAN FAILURE STRESS’ 

‘VECTOR, FSVM, AND TO THE APPLIED STRESS RATE VECTOR, ARV,’ 

‘ FOR DYNAMIC FATIGUE, FOR EXAMPLE' 

•TFVM— TFMl.TFMl,...’ 

‘ASV-AS\yAS2,...’ 

‘WHERE TFMl IS THE MEDIAN TIME TO FAILURE AT APPLIED ’ 

‘STRESS ASX, OR’ 

‘FSVM- FSM\,FSMl,...’ 

‘ARV-ARX^iR2,...’ 

‘ WHERE FSMl IS THE MEDIAN FAIL URE STRESS A T ’ 

‘APPLIED STRESS RATE ARV 


‘IF HOMOLOGOUS RATIO ANALYSIS IS BEING USED FOR ’ 

‘STATIC FATIGUE, THEN ENTER THE TIME TO FAILURE VALUES’ 




[44] 

[45] 

[46] 

[47] 

[48] 

[49] 

[50] 

[51] 

[52] 

[53] 

[54] 

[55] 

[56] 

[57] 

[58] 

[59] 

[60] 
[61] 
[62] 

[63] 

[64] 

[65] 

[ 66 ] 

[67] 

[ 68 ] 

[69] 

[70] 

[71] 

[72] 

[73] 

[74] 

[75] 

[76] 

[77] 

[78] 

[79] 

[80] 
[81] 
[82] 

[83] 

[84] 

[85] 

[ 86 ] 
[87] 


[ 88 ] 

[89] 

[90] 

[91] 

[92] 

[93] 

[94] 

[95] 

[96] 


‘ASA MATRIX TFX. THE ROWS OF TFX ARE THE VECTORS TFVH\,TFVH2,...ETC: 

‘ASSOCIATED WITH THE APPLIED STRESSES ASl y AS2,...ETC. AT WHICH THE ' 

‘FAILURE TIMES WERE MEASURED. THE ELEMENTS OF THE VECTORS TFVH(I )' 

'ARE RANKED IN THE ORDER OF INCREASING TIMES TO FAILURE. THE ’ 

‘DIMENSIONS OF THE VECTORS TFVH(I ) ARE THE NUMBER OF SPECIMENS TESTED ' 

‘AT EACH APPLIED STRESS. IF ALL THE SPECIMENS DID NOT FAIL,’ 

‘ZEROSM, ARE PUT IN FOR THEIR TIME TO FAILURE. IF THE' 

‘DIMENSION OF THE VECTORS TFVH(I ) ARE STILL NOT EQUAL DUE TO ' 

‘DIFFERENCES IN THE NUMBER OF SPECIMENS TESTED AT EACH APPLIED ’ 

‘STRESS, THEN MAKE THEM EQUAL BY FILLING THE LAST ELEMENTS OF THE ’ 

‘SMALLER DIMENSIONED VECTORS WITH ONES,(\). A SMALL DIMENSIONED ' 

‘ VECTOR MA Y THUS CONSIST OF A SERIES OF INCREASING TIMES TO ' 

‘FAILURE FOLLOWED BY A SERIES OF ZEROS FOLLOWED BY A SERIES OF ONES, FOR EXAMPLE-,’ 

‘TFVm — TFl,TF2,TF3,...0,0,0,...l,l,V 

•AS\ — VALUE OF THE FIRST APPLIED STRESS ’ 

‘ASV—AS\rAS2,...’ 

‘TFX— If PASV),( pTFVHl )) p TFVHl , TFVH2, ...’ 

‘WHERE TF1 IS THE SHORTEST TIME TO FAILURE AT APPLIED ’ 

‘STRESS /4S1 AND TF1 IS THE NEXT SHORTEST TIME ETC. OR,’ 

t > 

‘IF HOMOLOGOUS RATIO ANALYSIS IS BEING USED FOR DYNAMIC FATIGUE,' 

‘THEN ENTER THE FAILURE STRENGTH VALUES AS A MATRIX FSX.’ 

‘THE ROWS OF FSX ARE THE VECTORS FSVH1,FSVH2,.. .ETC.’ 

‘ ASSOCIATED WITH THE APPLIED STRESS RATE AR\^\R2,...ETC.' 

‘A T WHICH THE STRENGTHS WERE MEASURED. THE ELEMENTS OF THE ' 

' VECTORS FSVH(I ) ARE RANKED IN THE ORDER OF INCREASING ' 

‘STRENGTH. IF THE DIMENSIONS OF THE VECTORS FSVHfi ) ARE ’ 

‘NOT EQUAL, THEN MAKE THEM EQUAL BY FILLING THE LAST’ 

‘ELEMENTS OF THE SMALLER DIMENSIONED VECTORS WITH ZEROS, (0).’ 

‘FOR EXAMPLE-,’ 

‘FSVHl —FSl,FS2 FS22,0,0,0 ’ 

'ARI— VALUE OF THE FIRST APPLIED STRESS RATE ’ 

‘ARV—AR\AR2,...’ 

‘FSX—(( PARV),( p FSVHl)) pFSVHl,FSVH2,...' 

‘WHERE FS\ IS THE LOWEST FAILURE STRENGTH AT APPLIED’ 

•STRESS RATE ARI AND FS2 IS THE NEXT LO WEST FAIL URE ’ 

‘STRENGTH, ETC.’ 


‘AFTER THE DATA INPUT CONTINUE THE PROGRAM BY TYPING ’ 

‘THE COMMAND— STM,— DYM, — STH OR — DYH DEPENDING ON WHETHER STATIC’ 
‘FA TIGUE OR DYNAMIC FA TIGUE DA TA WITH A MEDIAN ANAL YSIS OR STA TIC ’ 
‘FATIGUE OR DYNAMIC FATIGUE WITH A HOMOLOGOUS RATIO ANALYSIS ARE ’ 
‘BEING USED TO GENERATE THE FATIGUE PARAMETERS N AND B ’ 

SIM— .((+A//)x« »2)+ « SO 


SIM = S. = e ^ s o + 75" ln ln 2 

i * 

I— 1 

7=1, control variable 
X2V— 10 
Y2V— 10 
XiV— 10 
Y3V— 10 

SIM: Y2 V— • TFVM 
Y2V = In t 
XIV— « ASV 
X2V = In o a 
SLS-.J3 — 2, pX2V 


initialize 

variables 

value 


slope intercept 

In tj. = -N In o fl + In B + {N- 2) In S. , text equation (7) with median values 
Y2V X2V 

In (r,S?) = -N In + In B , text equation (14) 

/ I no 

slope intercept 


73 = a two-element vector: first element 2, second element pX2V 


X2—J3 PX2V.Y2V 


X2 = a two-row, multicolumn matrix: first row X2V, second row Y2V 
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[97] STATIC*— 1 PCF XI 

Static = a two-element vector containing slope and intercept of text equation (7) or ( 14) 

[98] N— - STATIC[\\ 

N = —slope 

[99] —(I>l)/BHS 

Branch statement differentiating between calculation of B using homologous ratios (/> 1) or median values (/= 1). 

[100] B — *STA TIC[2] — (N—2)x « SIM 

B _ 2 ( intercept ) - (N-2) In S, 

[101] — (I=l)/VS 

Skip the next statement that calculates B for homologous ratio analysis and start calculating data variances. 

[102] BHS'.B — *STA 77C[2] 

B = e ( i ’" erce ’’b = e ln B 

[103] KS: VY2 — ( + /((Y2V- STA 7YC[2] +X2Vx STA 7'/C[l])*2)) 4- (( py2K)-2) 

J , 

VY2 = V(y) = (1/7—2) E \y. -(ax. + b)f 

i - 1 1 ‘ 

[104] SUM— + /X2V 

J 

SUM = E x 
l 1 

[105] VX2 — ( + /X2V-SUM+( pY2V))*2) + ( PYZV) 

VX2 = R(x) = (1/7) ,2, (Xi - x) 2 

[106] KSH — Ky2+KX2x( pKV) 

KS71 = K(a 2 ) or V(a 4 ) = 1/7 KO)//?(x) 

[107] VST2— VSTl x VX2 + (SUM+( PY2V))*2 

VST2 = V(b 2 ) or V(bJ = KSH [/?(*) -jE 2 ] 

[108] COVn— -VST\xSUM+( pY2V) 

COVI2 = Cov(a 2 ,b 2 ) or Cov(a 4 ,b 4 ) = -(KS71) x 

[109] — CT 

Go to instruction for formatting the graph function. 

[110] STH:R — PASV 

R = control variable equal to the dimension of the applied stress vector 

[111] SLP: TFVH— TFX[I\] 

TFVH = a vector whose elements are the Ith row of the TFX matrix 

[112] TFVH— (l * TFVH)/ TFVH 

TFVH is compressed to eliminate ones from the vector. 

[113] SVI—SV 

Define initial value of inert-stress vector SVI. 

[114] SL: — ((p TFVH) = ( p SVI ))/SE 

Branch statement comparing the dimension of the time-to-failure vector with that of the inert-strength vector. 

[115] — (( p TFVH) >( pSVI))/TFG 

Branch statement preparatory to making the dimension of TFVH and SVI equal. 

[116] DP—UpSVD 

DP = a random number between 1 and pSVI 

[117] SVI [DP] — 0 

SVI [DP] = 0 

[118] SVI— (0* SVI)/ SVI 

SVI is compressed to eliminate zeros. 

[119] — SL 

Go to the discrimination statement for the vector equilibration loop for static fatigue. 

[120] TFG-.DP - ?( p TFVH) 

DF = a random number between 1 and P TFVH 

[121] TFVH[DP] — 1 

TFVH[DP] = 1 

[122] TFVH— (l * TFVH)/ TFVH 

TFVH is compressed to eliminate ones. 

[123] —SL 

Go to the discrimination statement for the vector equilibration loop for static fatigue. 



[124] 

[125] 

[126] 

[127] 

[128] 

[129] 

[130] 

[131] 

[132] 

[133] 

[134] 

[135] 

[136] 

[137] 

[138] 

[139] 

[140] 

[141] 

[142] 

[143] 

[144] 

[145] 

[146] 

[147] 

[148] 

[149] 

[150] 


SE: TFVH—( 0 * TFVH)/TFVH 

TFVH is compressed to eliminate zeros. 

SVI—( pTFVH)\SV 

SVI is shortened to the first pTFVH elements. 

Y2V—Y2V,® TFVHxSVDl 
y = In (t f S 2 ) 

X2V— X2V,% ASV[I]+SVI 
x = In a Q /S j 

/— 7+1 

Increase control variable by one. 

— (I<R)/SLP 

Branch statement to go through loop again ( I<R ) or go to least squares routine (/ >R) 

— SLS 

Go to least squares routine for static fatigue. 

DYM:Y2V— « FSVM 
YiV = y = Ini 
XiV— ® ARV 
Xi V = x = In a 
DLS:Ji — 2, pXiV 

J3 = a two-element vector: the first element of which is two and the second element pXiV 
Xi—Ji PX3 V,Y3V 

Xi = a two-row, multicolumn matrix: the first row X3V, second row YiV 
D YNAMIC — 1 PCF Xi 

DYNAMIC = a two-element vector containing the slope and intercept of text equation (5) or (13) 
W~(+DyAMM/C[l])- 1 
N = (1/slope — 1) 

-~(I>1)/BHD 

Branch statement differentiating between the calculation of B using homologous ratios (/> 1) or median values (/= 1) 
B'~-*((N+\)x.DYNAMIC[2]) — (to N+l) + (N—2)x ® SIM 

B = e (W+l )y.(intercept) - [In (N+l) + (N- 2) In 

— (I=1)/VD ' 

Skip the next statement that calculates B for homologous ratio data and start calculating data variances 
BHD:B—*((N + \)xDYNAMIC[2\) - « N + 1 

B = g(W+l) x (intercept) - In (N+ 1) 

VD: VY3 ■— ( + /((Yi V- DYNAMIC^} +Xi VxDYNAMIC(\])*2)) + (( pY3V)-2) 

J 

VYi = (l/J-2) £ \y -(ax.+b)] 1 

i= l ' ' 

SUM" — V/XiV 
] 

SUM = E Xj 

i= 1 

VXi~(+/(XiV-SUM+( pYiV))*2)-i-(pYiV) 

J 

VXi = R(x) = (1/7) E (x, - x) 2 

i = i 

VDYl «— VY3 + VXi x( pYiV) 

VDY1 = K( 0] ) or V(a } ) = (1/7) V(y)/R(x) 

VDYl— VDYl x VXi + (SUM +( py310)*2 
VDYl = K(6,) or V(bJ = VDYl(R(x)+f 1 } 

COV12—-VDYlxSUM+( pYiV) 

COV1 2 = Cov (ffpii,) = Cov(a 3 ,b } ) = ~(VDYl)x 

— CT 

Go to instructions for formatting the graph function. 

DYH.R — pARV 

R = control variable whose value equals the dimension of the applied stress vector 
DLP:FSVH—FSX[I,} 

FSVH = a vector whose elements are the Ith row of the FSX matrix 
SVI—SV 

SVI = SV defines initial value of the inert-strength vector 
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[151] FS VH— (0 ± FS VH)/FSVH 

FSVH is compressed to eliminate zeros. 

[152] DL: — * (( p FSVH) = ( P SVI))/DE 

Branch statement comparing the dimension of the failure stress vector with that of the inert-strength vector. 

[153] — (( pFSVH)>{ pSVWFSG 

Branch statement preparatory to making the dimension of FSVH and SVI equal. 

[154] DP-l(pSVI) 

DP = a random number between 1 and pSVI 

[155] SF7[DP]~ 0 

SVI[DP] = 0 

[156] SVI— (0* SVI)/ SVI 

SVI is compressed to eliminate zeros. 

[157] 

Go to the discrimination statement for the vector equilibration loop for the dynamic fatigue. 

[158] FSG:DP—T( p FSVH) 

DP = a random number between 1 and pFSVH 

[159] FSVH[DP] — 1 

FSVH[DP] = 1 

[160] FSVH— (1 * FSVH)/ FSVH 

FSVH is compressed to eliminate ones. 

[161] — DL 

Go to the discrimination statement for the vector equilibration loop. 

[162] DE:nV— Y3V, * FSVH + SVI 

Y3V = y = In S/S f 

[163] X3V—X3V, « ARVII]+SVM 

X3V = x= In b/Sf 

[164] 7—7+1 

Loop control variable update. 

[165] —(I<R)/DLP 

Control variable comparison branch statement. 

[166] —DLS 

Go to least squares routine for dynamic fatigue. 

[167] CT:‘ DECIDE ON THE LENGTH OF THE LN( APPLIED STRESS) AXIS BY ’ 

[168] ‘ASSIGNING A VALUE TO THE VECTOR LASV. FOR EXAMPLE ’ 

[169] ‘LASV— 16.45 + 0.05 x 140’ 

[170] ‘BEGINS THE AXIS AT 16.50 WITH A SCALE OF 0.05 FOR 40 ’ 

[171] ‘POINTS, THUS ENDING THE AXIS AT 18.45 ’ 

[172] 

[173] ‘NEXT DECIDE WHICH TYPE OF DESIGN DIAGRAM IS DESIRED.’ 

[174] ‘FOR MULTIPLE PROBABILITIES WITHOUT CONFIDENCE LIMITS,’ 

[175] ‘SPECIFY FAILURE PROBABILITIES AND PERFORM THE FUNCTION ’ 

[176] ‘LTF AFTER EACH SPECIFICATION, STORING THE RESULTIN' 

[177] ‘AN APPROPRIA TEL Y NAMED VARIABLE, FOR EXAMPLE ’ 

[178] ‘F*-.oor 

[179] ‘LTFOOl-LTF’ 

[180] 'FOR A DESIGN DIAGRAM AT A SPECIFIC FAILURE PROBABILITY’ 

[181] ‘ WITH CONFIDENCE LIMITS SPECIFY THE FAILURE PROBABILITY ’ 

[182] 'AND PERFORM EITHER THE FUNCTIONS LTF AND SVLTF OR LTF ’ 

[183] ‘AND DVLTF AFTER THE SPECIFICATION, STORING THE RESULT IN AN’ 

[184] ‘APPROPRIA TEL Y NAMED VECTOR, FOR EXAMPLE ’ 

[185] ‘F— .001 ’ 

[186] ‘LTFOOl-LTF’ 

[187] ‘DVLTF001 —DVLTF’ 

[188] 

[189] ‘NEXT ENTER THE ARGUMENT FOR THE GRAPH FUNCTION, FOR EXAMPLE ’ 

[190] ‘X—l 40 PLASV,LTF00l,LTF0l,LTFl,LTFS,LTF99,LTF999 ’ 

[191] ‘OR’ 

[192] ‘C90*-1.645xDKL77 : 001*.5’ 

[193] ‘LL—LTF001—C90’ 



[194] 

[195] 

[196] 

[197] 

[198] 

[199] 

[ 200 ] 
[ 201 ] 
[ 202 ] 
[203] 


[ 1 ] 

[ 2 ] 

[3] 

[4] 

[5] 

[ 6 ] 

[7] 

[ 8 ] 

[9] 

[ 10 ] 
[111 
[ 12 ] 


[13] 


[14] 


[15] 


[16] 


[17] 

[18] 


[19] 


[ 20 ] 

[ 21 ] 

[ 22 ] 


‘ UL — LTF001 + C90’ 

‘X—4 40 pLASV,LL,LTF001 , UL' 


•IN ORDER TO GENERATE A TABLE OF GRAPHED DATA USE ' 

•THE TRANSPOSITION FUNCTION I?, FOR EXAMPLE ' 
'V—LASV,LTFOOl,LTPOl,LTFl,LTF5,LTF99,LTF999 ’ 

'Kl — 7 40 pV' 

‘$K1 ’ 

•FINALLY TO CONTINUE THE PROGRAM TYPE THE COMMAND — Cl ' 

C2.GRAPHX 

V 

VSVLTF\U\V 

WLTF—SVL TF;KI;A 1 ;A 2; A 3;KS;LA S VI ;A LA S;RLA S; VN;CO VNLB; VLB; A 1;A &;A 9; A 1 0;S/;71 ;72;70 
‘IF YOU ARE USING MEDIAN ANALYSIS AND ASSUMING WEIBULL’ 

' DISTRIBUTIONS FOR LOG TIME TO FAILURE AND INERT STRENGTH' 

•CONTINUE THE CALCULATION BY TYPING THE COMMAND — W.' 

'IF YOU ARE NOT ASSUMING WEIBULL DISTRIBUTIONS AND ARE USING ’ 

•ONL Y FA TIGUE AND STRENGTH DA TA TO DEFINE THE ERROR LIMITS ON ' 

•fatigue'* PARAMETERS AND LOG INERT STRENGTH, THEN CONTINUE THE ' 

• CALCULATION BY TYPING THE COMMAND —M FOR MEDIAN ANALYSIS OR ' 

• — H FOR HOMOLOGOUS RA TIO ANAL YSIS ' 

W:J 1 —pSV 

J { = dimension of the inert-strength vector 
J2 — PASV 

J 2 = dimension of the applied stress vector 
70*- 71 x72 

J o = y i x '2 

KI— ((MI*2) x ((\2 + MI) -s- (1 1 h- A/ 7)*2)-l) + ((ft ft 2)-MIx ft !l-=-A/7)*2 
K. = m f (1 +2/m ? ~ r2(1 + [i n l n 2 -m In T(1 + 1/m)] 2 

' ' r 2 (l + l/mp ' f 

.41— ((ft ft 2)-(MI+N-2)x ft !2 + (A7-2)+MQ*2 


a 3 = r Y 1 + ~~ 
m. 


A, = [in 1“ 2 - ln T(1 +— )] 

L N-2 m,- J 

A2 *- !2 X (N-2) +MI 

, 2 ( 77 - 2 ) . 

.4 3 — (1 (N-2) + MI)*2 

X 

KS—(((MI+N-2)*2) X (A2+A3)-l) +A 1 

* ’ <"/«’’[ r ( 1 '^) - ! ( -f )]/ ** -f )♦ [ ■” 

LASVl — ft ASV ' ' ' 

LASVl = ln o a 
ALAS—( + J2)x + /LASVl 
1 h 

ALAS = ln a Q = — E ln a a j 

RLAS— ( + 72) x + /(LASV1-ALAS)*2 

1 J 2 2 

RLAS = 7?(ln a a j) = — E (In a a j - ln afy 


J 2 i= 1 


VN— (JKS X (N- 2)*2) - 70 X (MI*2) x RLAS 
VN = V(N) = K s (N-2) 2 /J 0 m. 2 7?(ln <rj 
CO VNLB— VNxALAS 


COVNLB = Cov(N, ln B) = V(N) IrTa 
VLB— (VNx ((ALAS- ft SIM)*2) + RLAS) + ((N-2)*2)xKI+J\ xMI*2 
VLB = V (ln B) = V(N)((ln - InT) 2 + fl(ln aj( + (N-2) K/J ] m 2 
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[ 23 ] 


Al — (U2+MI) + (ll-i-MI)*2)-l 

A 7 = [r(l +2/m j ) - r 2 (l + l/m.)]/r 2 (l + 1/m.) 

[24] AS-'— ® « 1-hI-F 

A g = In In (l/(l-7)) = x p 

[25] A9 — m \l+MI 

A 9 = In T(1 + 1/mp 

[26] A10—(A8*2)-v-Jl xMI*2 

A io = x F 2 ( 1 / 7 , m?) 

[27] VLSI— A10 + ((Al + A9*2) + Jl)-2xA8xA9 + JlxMI 

VLSI = K(ln S.) = x 2 (l/7j m. 2 ) + 1/7, [{[r(l + 2/m ; ) - r 2 (l + l/mp]/r 2 (l + 1/m.)} 

+ [lnr(l + 1/m.)] 2 ] - (x p 2/7, m.) In T(1 + 1/m.) 

[28] — 'C3 

[29] M:VN—VST\ 

VN = V(N) = KS77 = V(a 2 ) 

[30] VLB ■'—((( ® SIM)*2 ) x KS71) + VST2 + (2 x ( ® S/A7) x COV\2) + ((2+STATIC[\})*2 ) x VLSI^ ( p S10 

VLB = K(ln B) = (In 3.) 2 K(a.) + K/6,) + 2(ln S.) Cov(a,A) + ((a, + 2) 1/7.) K(ln S.) 

[31] COVNLB — -COV12 

COVNLB = Cov(N, In B) = -Cov{a 2 ,bJ 

[32] — C3 

[33] H:VN—VST\ 

V(N) = VSTl = K(a 4 ) 

[34] VLB— VST2' 

VLB = F(ln B) = KS72 = V(bJ 

[35] COVNLB— -CO VI 2 

COVNLB = Cov(N, In B) = -Cov 12 = -Covfa^ bj 

[36] C3:S/— *((-*- AW) x ® ®(1 -i- 1-7)) + ® SO 

Si = s = e In ln + In S 0 

[37] AS — *LA,SK 

[38] 'VN EQUALS’ 

[39] VN 

[40] ‘ VLB EQUALS’ 

[41] VLB 

[42] * COVNLB EQUALS’ 

[43] COVNLB 

[44] ‘VLSI EQUALS’ 

[45] VLSI 

[46] VLTF— VLB + (VN X ( ® SI+AS)*2) + ( VLS!x(N-2)*2) + 2 X CO VNLB X ( & SI+AS) 

VLTF = K (ln iy) = K(ln B) + (ln S/af V(N) + (N-2) 2 K(ln S.) + 2 Cov (N, ln B) ln (SJoJ 

VDVLTF\U]V 

VVLTF—DVLTF;AhA2;A3;A4;KD;LARV;ALAR;RLAR;A5;VN-,A6;CONVNLB;VLB;Al;A8;A9;A10;JUJ2;J0;SI 

[1] ‘IF YOU ARE USING MEDIAN ANAL YSIS AND ASSUMING WEIBULL ’ 

[2] ‘DISTRIBUTIONS FOR LOG TIME TO FAILURE AND INERT STRENGTH, ’ 

[3] ‘CONTINUE THE CALCULA TION BY TYPING THE COMMAND — W. ’ 

[4] 'IF YOU ARE NOT ASSUMING WEIBULL DISTRIBUTIONS AND ARE USING ’ 

[5] ‘ONL Y FA TIGUE AND STRENGTH DA TA TO DEFINE THE ERROR LIMITS ON ’ 

[6] ‘FATIGUE PARAMETERS AND LOG INERT STRENGTH, THEN CONTINUE THE ’ 

[7] 'CALCULA TION B Y TYPING THE COMMAND — M FOR MEDIAN ANAL YSIS OR —H FOR ’ 

[8] 'HOMOLOGOUS RA TIO ANAL YSIS. ’ 

[9] W:J\— pSV 

7, = dimension of the inert-strength vector 

[10] J2— pARV 

J 2 = dimension of the applied stress vector 

[11] JO— 71x72 

J Q = 7j x 7 2 = total number of specimens 

[12] Al — l+(N-2)-i-MIxN+l 

A t = 1 + (N-2)/m. (N+ 1) 
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[13] >12 — 1+2x(N-2)xMIxN+1 

A 2 = 1 + 2 (N-2)/m.(N+ 1) 

[14] A3~-MIx(N+l)+N-2 l 

A 3 = m. (N+ \)/(N-2) 

[15] >14 — (( & & 2)-A3 x « Ml-1)*2 

/1 4 = {In In 2 - m. [(//+ l)/(AT-2)] In T [1 +(N-2)/m. (N+ 1)]} 2 

[16] KD— >14 + 043.2) x ((042-1) + (!>11-1)*2)-1 

K d = m ) I(N+ 1 )/(W-l)] 2 {r (1 + (2/mfiN-2)/(N+ 1)) - r 2 (1 + (N-2)/mfN+ l))}/r 2 (l + (N-2)/m.(N+ 1)) 

[17] LARV— ® ARV 

LARV = In b 

[18] ALAR — (.+J2)x +/LARV 

h 

ALAR = In a = (\/J.) E In b, 

2 j= i y 

[19] RLAR — ( + J2)x + /(LARV-ALAR)*2 

J 2 _ , 

RLAR = R (In a) = (1/7.) E (In a,- In a) 2 
z y-i y 

[20] t45 — ((N-2) X (N+ 1))*2 

/lj = (N-2) 2 (N+ 1) 2 

[21] VN—KDxA5 + J0xRLARxMl*2 

VN = V (N) = K d (N-2) 2 (N+ l) 2 /7 Q m 2 i? (In a ) 

[22] >16 — ( d BxN+\)+ALAR-\ + 3x © S/M 

= In [B (W+ 1)] - 3 In l - 1 + TiT* 

[23] COVNLB—(VN+N+ l)xALAR 

COVNLB = Cov (N, In B) = KJ(N- 2) 2 (N+ 1) In a/J Q m 2 . R (In b ) 

[24] VLB—(VN+ (N+ 1). 2) x (RLAR x 1 + 72) + (-A6)*2 

VLB = V (In B) = {K d (N-2) 2 /J Q m 2 R (In b )} {[3 In 5 . + 1 - In (B (N+ 1)] - Tn & ] 2 + (1 + J 2 ) R On 6)} 

[25] /17-*— ((!2-i-A/7)-t-(!l -s-A7/)*2) -1 

A 7 = [r (1+2/m.) - r 2 (l + l/m.)]/r 2 (l + 1/m.) 

[26] >48 — ® ® Ih-1-F 

A s = In In [1/(1— 7*)] = x p 

[27] >19 — & ll+Mf 

A g = In T (1 + 1/mp 

[28] >110— (.48*2) -i- 71 xA//»2 

' A 10 = 4 V{ - Um ) = X F /J 1 4 

[29] VLSI—A10 + ((A7 + A9*2) -s- 71)-2 x>18 X>19-r71 x MI 

VLSI = V (In S) = x 2 V (1/mp + V (In S Q ) + 2 Cov (m. , In SJ x p 

[30] — C3 

[31] M-.VN— ((+DYNAMIC[\])*4)xVDY\ 

VN = V (N) = (l/a,) 4 K (op 

[32] >1 12 — (-(Z?y//>1A//C[2]) -s- (DYNAMIC[1])*2) + (+ DYNAMIC[l]) + ((-i-DYNAMIC[\])*2) x & S/A/ 

/ll 2 = (-6,/a 2 + 1/a, + On §p/a 2 ) 

[33] VLB — (((/4 1 2)*2) X KOri)+(((-ByM4A//C[l])*2) x K/)J7) + (2 x 0412) x (^DYNAMIC [\]) x COK12) + (((^ DYNAMIC[\])-3)*2) 
xVLSI+(pSV) 

VLB = V (In B) = (-b { /a 2 + 1/a, + (1/a 2 ) In £.) 2 K(a,) + (1/a,) 2 K(6,) + 2 (- 6,/a, 2 + 1/a + (1/a 2 ) In Sp(l/a ) 

Cov (a, ,6, )+ 1/7, (1/a, - 3) 2 K (In Sp ' . i 

[34] CO VNLB — ~(N+ 1) X ((+DYNAMIC(\])*2) x COVi 2 

COVNLB = Cov (N, In E) = - ( N+ l)(l/a 2 ) Cov (a, ,6,) 

[35] — C3 

[36] //: VN— ((-*-/) YNAMIC{ 1 ])*4) x K/m 

0W = K(N) = (l/aj 4 ) V (a 3 ) 

[37] >113 — (-(D m4 A//C[2]) -h- (/>yAMA//C[l])*2) + (^Z)y//>lA//C[l]) 

>1 13 = (~b } /a 2 + l/a 3 ) 

[38] KiB- (((-413)*2) X VDY\) + (((+DYNAMIC[\\*2) x VDY2) + 2 X 04 13) X (+ DYNAMIC[\]) x COVY2 

VLB = V (In B) = (- b/a^ + 1 /a } ) 2 V(a 3 ) + (1 /af V (b 3 ) + 2 (-ya 2 + l/a 3 )(l/a 3 ) Cov (a y b 3 ) 

[39] COVNLB— ~((+DYNAMICll])t2)xCOV\2 

COVNLB = Cov (N, In B) = (-1/a 2 ) Cov (a 3 ,b 3 ) 
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[40] C3:S/— *(+MI)x ® ® (l-I-i=))+ ® SO 

- (l/m,) In In (1/1-f) + In S 

S. = e 0 

[41] /IS— *LASV 

[42] 'VN EQUALS’ 

[43] KIV 

[44] * VLB EQUALS' 

[45] VLB 

[46] ' COVNLB EQUALS’ 

[47] COVNLB 

[48] 'VLSI EQUALS' 

[49] VLSI 

[50] VLTF— VLB + ( VN X ( 6 SI+AS)*2) +(VLSIx (N-2)*2) + 2 x COVNLB X ( ® S/-/1S) 

KLTF = V (In t f ) = V (In B) + On S/of V (N) + (N-2) 1 V On S.) + 2 Cov (N, In B) In (S/aJ 

VLTFl D]V 
VR—LTF;SI 

[1] S/— *(( + MI)x ® 6 (1 + 1-F)) + ® SO 

(l/m:) In In (1/1-f) + In S„ 

S/ = S ( = e 0 

[2] 7J — (-7Vx/^4S10 + ( ® B) + (7V-2)X ® S7 

In /, = - TV In a + In B + (N-2) In S. 

/ a t 

V 


Example of the Use of the Program 

The program DESIGN was used to initiate the iterative process of hardware development for the meter square 
spark-chamber frames for the EGRET. The frames were made from Comings 9658 machinable glass-ceramic, 
Macor™. Reports on the static and dynamic fatigue testing of this material appear in references 9 and 10, respec- 
tively. This section of the report contains design diagrams for Macor, using a homologous ratio analysis of dynamic 
fatigue testing data. Two diagrams are presented: One diagram (Figure 1) is a plot of equation (11) for probabilities 
of failure of 0.001, 0.01, 0.1, 0.5, 0.99, and 0.999. The other (Figure 2) is a plot of equation (11) with error bars 
at the 90-percent confidence level for a probability of failure of 0.001. 

The homologous ratio analysis of static fatigue data for this material was interesting. A plot of the homologous 
ratio data is shown in Figure 3. The figure indicates that there is an anomaly in a graph of the data around In afS. 
equal to -0.52. The anomaly is the rapid increase in In {tj S 2 ) at values of In of S. below -0.52. For a material that 
is susceptible to subcritical crack growth, such a plot is expected to be a straight line with a slope whose negative 
is the fatigue parameter N (equation 14). The interpretation of this anomaly is that it represents a fatigue limit in 
this material. 


17 



LN APPLIED STRESS IN PASCAL 










52.00 





or 

45.00 




38.00 


-0.600 


1 

-0.500 


I 

-0.400 
In {<7 a /Sj) 


-0.300 


- 0.200 


Figure 3. Static fatigue homologous ratio least squares plot for Macor. 
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APPENDIX A 


SOME EFFECTS OF A DISTRIBUTION IN THE SIZE OF 
THE LARGEST FLAW IN A MATERIAL 


The strength of high modulus, low K JC materials is very sensitive to surface damage. A plot of failure strength, 
a , versus the square root of the failure-initiating flaw size, a'f, is a hyperbola with the axes as the asymptotes. A 
plot for K JC equal to values ranging from 0. 1 to 10 MPa m‘ A is shown in Figure A-l . The plot is intended to illustrate 
that proof-testing is normally necessary to ensure the quality of glasses and ceramics and that improving surface quali- 
ty to get higher strengths does not necessarily improve reliability as discussed below. 

Normal nondestructive test techniques do not consistently find flaws less than 1 mm (a' A < 0.03 m' /j ). For metals 
(K IC » 10 MPa m 1/2 ), this is normally no problem because corresponding strengths are much greater than 150 MPa 
and a 10-percent variation in the size of the worst flaw in a piece of material does not appreciably change observed 
strength. However, for a glass or ceramic material (K /c ~0.5 to 5 MPa m /2 ), flaws > 0.1 mm (a^ 2 > 0.01 m /j ) are 
at the lower size limit of even special NDE, and such flaws result in strengths which are unacceptably low. In addition, 
a 10-percent variation in the size of the worst flaw in a piece of such material has a major effect on the observed 
strength of that material since observed strengths lie on the high strength leg of the K JC hyperbola. Therefore, im- 
proving processing parameters to increase average strength does not necessarily improve reliability. 
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Figure A-l. Failure stress versus square root of crack length for critical stress-intensity factors 

( K IC = 0.1, 0.5, 1, 5, and 10 MPa m' /j ). 
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APPENDIX B 


DERIVATION OF THE DYNAMIC FATIGUE EQUATION 

(after reference 19) 


V = -^-= AK ", but K, - Yo\fa\ therefore, V = A(Ya\Ja)" 
at 1 1 

Rearranging factors, da/a" n = A Y"o" dt and recognizing that b = do/dt is a constant and that dt = da/ a results 
in the integral expression 


s 


V 


da 

a"' 2 


\° f AY" o" do 

a. a 

I 


which on evaluating gives 


-(<¥) 


2a 


(2 -AO 


a / 


AY"o 


N + 1 


b(N+ 1) 


\ a f 


Iff; 


/V-2 _ N-2 

Neglecting a' - " 2 " relative to aT ~ T ~ and assuming o. = 0 gives 


2 -(-^) ziy^ff^ 1 

(N-2) a '- 2 = (N+ l)j 

But a. is the starting flaw size in the absence of slow crack growth. It is related to inert strength S. by the expression 
a ; = ‘k 2 ic /S 2 . Y 2 . Thus 


2 , AY N o N f + 1 

(N-2) \ s 2 y 2 / = (N+1) or 

i 


but Oj. is observed strength. 
Therefore, 


where 


S" +l 


2 (N+ l)ff S?- 2 
(N-2) A Y 2 K^ c 2 


B (N+ 1) a Sf 2 


B = 2/ (N-2) A Y 2 K"~ 2 

This equation can also be used to estimate how fast a stressing rate, b, must be used before subcritical crack growth 
is eliminated and inert strengths are being measured. The inert-strength stressing rate is 

a. = S 1 ./B(N+ 1) 
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APPENDIX C 


DERIVATION OF THE STATIC FATIGUE EQUATION 

(after reference 23) 


V = — = AK " but K, = Ya Va and a = a , which is a constant. Thus, 
dt 1 1 a 

V = A(Ya Va f or = AY N ^dt 


which on integrating gives 


j V ■%= 


N- 2 


N-2 


Assuming t. = 0 and neglecting a^ 2 relative to a. 2 gives, on evaluating the integrals, 

2a^ 2 

(N-2) = AY* o%t f 


but again as with dynamic fatigue, 


Therefore, 


Solving for tj gives 


a, = K)c/S? Y> 


N-2 


K\ 


(i^)(F#) 


__ 2 S?-\- N 

~ AY\N- 2) K”- 2 


= BS?- 2 o- n 


B = 


AY\N- 2) ^- 2 




APPENDIX D 


ADDITIONAL CONCERNS IN DESIGN-ALLOWABLE DEVELOPMENT 


For materials that do not possess a characteristic strength such as is commonly associated with metals, developing 
design allowables for a piece of hardware is an evolutionary process between design and materials engineers. It typical- 
ly starts with an initial design with a certain stress distribution containing a maximum stress that is checked relative 
to the applied stress given by equation (11) for the required mission lifetime. As it stands, however, equation (11) 
does not take many things into account that, when included, typically reduce the maximum allowable design stress. 
It assumes that the hardware will be a relatively small 3- or 4-point bend specimen. To correctly apply equation (1 1) 
requires that the anticipated parent flaw population that will cause failure in the hardware is the same as the test specimens 
(i.e., failure will occur from similar populations of edge, surface, or volume flaws). In addition, differences in the 
stress distributions between the test specimens and the hardware must be accounted for. Finally, the size difference 
between the two systems must be accounted for. An example of how this was accomplished for EGRET spark-chamber 
frames is given in references 9 and 10. Other useful information in this area is found in references 24 through 26. 

As previously noted, this type of analysis does not obviate the need for proof-testing. Proof-testing establishes 
a starting point for lifetime estimates based on crack dynamics derived from fatigue testing. Given proof-test levels 
and crack propagation data, margins of safety can be calculated. 
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APPENDIX E 


DERIVATION OF THE VARIANCE OF In S, 

(after reference 3) 


In S. = — x P + In S 
‘ m. F ° 

l 

x F : In In ( 

F = the cumulative failure probability = 

/ = rank from weakest to strongest 
n = number of specimens tested 


Applying the error-propagation rule to determine the variance in In S. in terms of the variability of m. and S o and 
assuming that x p is a constant gives 


,3 In S. 
V (In S ; ) = ( t - 


<d 1/m. 


\ 2 ,3 1ns. , 3 In S. w 3 In S„ 

) V *' m > + (jlTsO + 2 Cov( 1 / m ,,lnS,)( Fr ^)( 


3 In S 


}) 


= ** V{\/m) + V (In S o ) + 2 Cov (1/m,., In S Q ) x f 
but from reference 3, page 804, equations (A15), (A26), and (A27), respectively, 

V (1/m.) = \/J x m] 


1 r r(l + 2/m.) - r 2 (l + l/m.) , , i 

’'any =T, { r^l + Wmf I'nro + '^l 2 } 


1 

Cov (1/m., In S) rCov (m., In S ) 

1 0 m. to 


l 


1 m. 

7 In T (1 + 1/m,.) 


In the derivation of the foregoing three expressions, the authors in reference 3 assume that x p is an independent 
variable, making equations (20) through (22) applicable to the evaluation of V (m,.) and V (In S q ). For what the authors 
accomplished in reference 3, this is a reasonable assumption because they generated fatigue data for an “ideal” material 
on a computer using a Monte Carlo technique. It was assumed that the fatigue parameters N and B were known and 
that the inert-strength distribution of this ideal material was given by a two-parameter Weibull distribution, the modulus 
and scale parameter of which were known. With the Monte Carlo technique, a given number of specimen strengths 
were then calculated by randomly selecting their failure probability from a uniform distribution between 0 and 1, 
making x p the independent variable. The result is that V (In S Q ) and Cov (m. In S o ) are independent of S and S o , 
depending only on the number of specimens J x and the Weibull modulus m,.. 

It can be argued that in the actual determination of m. and S Q , x f is not an independent variable since it is arrived 
at by simply ranking the measured strength values. The quantity V (x f ) = (lA/j-1) E(x f - x^ 2 could as easily belong 
to a sample of strengths with a range from 10 to 20 MPa as from 20 to 30 MPa. If measured strength values are 
considered, the independent variable and the linear expression of the Weibull equation are rearranged to reflect this. 
The result is 

Xp = m i In S,. - m,. In S g 
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Subsequently, going through an analysis similar to that in reference 3 gives expressions that differ from the ones derived 
previously. Because of this discrepancy and the desire to stay as close to the actual data as possible, error-bar estima- 
tion without assuming Weibull distributions is preferred. 

A technique sometimes used to do homologous ratio analyses with fatigue samples in which the number of specimens 
does not match the inert-strength sample size is to calculate Weibull parameters and then generate strengths correspond- 
ing to the rankings in the largest sample size. Since this method generates more data points than have actually been 
measured, it is not preferred. Using only the test data measured to determine the variances of the fatigue parameters 
without assuming any form of distributions is preferred. 
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APPENDIX F 


DERIVATION OF COVARIANCE EXPRESSIONS 


Differential Expression 

Reference 3, page 805, indicates 


Cov (N, In B) = V (TV) 


This expression can be derived from the definition of covariance in Bevington (reference 20, page 59) as follows: 

Given this definition, it is possible to arrive at a differential expression for covariance as follows: 

= lim — E [(^Wp) 2 ] = * im ±E (v.-v) 2 
/*v 7— oo J L\ v ' J J— oo J dv t ' 

= t, lim \L d lL)( “ m 4-E(^) ! ) = — VM 

V/— 0° ./ cfr./Vj—oo 7 ' / dv v ^ ’ 


Linear Regression Expression 

Reference 3, page 803, indicates 


Cov (cr,b) = - 


a 2 _ 


JR (x) 


This expression is derived similarly to the differential expression for covariance: 

lim 1„ \( b r b \,~ "\ 2 1 I™ A.^ db i f „_^2 


■i - [(£) 7 s 


But from equation (19), 


thus, 


and 

But 

thus, 


b = -ax + y 
db. 

I _ 

da. ' 
lim 1 


2 -lim „ . -, 2 _ / - hm 1 „ \/ lim 1 , _. 2 \ 

a 2 . = r —Ex {a. -ay ~ l T ~rE x. If T „ — (a -a) I _ 

ab 7— oo 7 i K i ' \ 7— oo 7 ' / V7— co J i / 

o 2 


- x r (a) 


V(a) = 


JR(x) 


Cov ( a,b ) = a 2 . = 


°b JR (. X ) 
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APPENDIX G 


DERIVATION OF DYNAMIC AND STATIC FATIGUE MEDIAN VALUE 
ERROR LIMITS NOT ASSUMING A PARTICULAR DISTRIBUTION 


Dynamic Fatigue 


In 


£ = ( ) In a + -) [In B + In (N+ 1) + (N-2) In £] 

y = a x x + b x 

‘,-jvT! • 4 1 = (^ T )HnB + ln(N + l) + (W-2)lnJ,] 

y = In % , x = In a 

given V (a.), F (6,, and Cov (a,, 6,) from equations (20) through (22). 

Solving for N in the equation for a, gives 
N = — - 1 

Applying the error-propagation rule, 

= (t; / = (i;) 4 ^".' 

Solving for B in the equation for b l gives 

In B = (N+ 1) b t - In (. N+ 1) - (N-2) In S. = + In a, - (i- - 3 ) In §. 


Applying the error-propagation rule, 

3 In B 2 


yon* - (^ )■ ^,> ♦ W ^ + * (*Sf K^f) 


d In B \ / d In B • 


+ titi)yon*> 

=(-l + ^ + s; ftof <y ,, <“.> + (-^) j ''».) 
+ 2 ( ■ % * i + ln s ‘ )(^) Cm ia '’ b,) 

+ (i--3)V(l„S,) 
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Strictly speaking, it is not possible to evaluate V (In S.) without either assuming a particular distribution or going 
through a parameter analysis such as the one performed for the Weibull distribution in reference 3 or else testing 
many samples of specimens. However, since V (In S.) is probably close to V (In S) for most distributions and since 
V (In S) = (\ //,) V (In S.) (Bevington, reference 20, page 21), the equivalence ofV(ln S.) and K(ln S) is assumed 
so that a particular distribution does not need to be specified when using only measured data to define V (In B). Thus, 


F,l„ B) = ( - ^ + 1 + -L In S, ) V (a,) + (-1) ! y (bi) 

+ 2 ('^ + ^ + ^ lnS < )(^) C ° V W 


Applying the differential expression for covariance, 


Cov (N, In B) 

d In B 
dN 

d In B 
db y 

d In B 
d In §. 

I 


Cov ( N , In B) 


d In B 
dN 


d In B dN 8 In B db, a In B 3 In $. 


( o m n aiy 

— 

dN dN 


+ 


dN dN db, dN 3 In S. dN 


■) V (N) 


(».-wft -■■*.) • 10= 1 

N+i • JZ’M -7CT(*.-NTi- |n ^) < NweI > 


a In s. 


db 


(N+ 1) — i V (A) = (N+ 1) Cov (b.,N) 
dN 1 


however, 


Cov = (6,,N) 


1 

- ~2 Cov (<*,»&,) (Note 2) 


Thus, 


Cov (N, In B) = - (N- 1) (-jj ) Cov (a,, 6,) 

Note 1: 


d A 

dN 


db i j 1 \f u 

dN ~ VAmT/v bl 


N+ 1 


■- In 



is obtained from the expression for the total derivative of b { = b { (N, In B, In 5.) 


36 



db x 

In 


36, dN 36, 3 In 5 3 In 5 3 In S. 

dN m + 3 In B dN 3 In 8. dN 


36 , / 1 \ / 1 \ 

— - + | )( 6. In S. ) 

3 N \N+ 1 /V 1 N + 1 '/ 


as 


3 In S. 

'■ = 0 

dN 


This differs from Ritter et al. (reference 3). On page 804 of that reference, the authors indicate that 36/3N = 
db /dN, which is true for an intercept that depends only on N, but in this case there is a dependence^ In B and 
In S, that, strictly speaking, should not be left out. The result is that the terms In B (N+ 1) and -3 In S. in equation 
(A43) of that reference drop out of the expression for covariance when a Weibull distribution is assumed as shown 
in Appendix H. 


Note 2: 


Cov (6,,iV) = -l/a\Cov (o,, 6,) because applying the error-propagation rule to^the expression for In B as a func- 
tion of o ,6,, and In S.; and applying it again for In B as a function of N, 6,, and In S.; and comparing terms show that 


* (^)(^) = 2 (> V ? ^ Cov ( °' A> 


must equal 


2 (^)(^r) C ° V = (V N> ■= ^ (N+ 1) ( 6, - ^ - in S, ) Cov (»,. 

= 2 (i-) (6, - 0 ,-/ 1 ! S.) Cov (6,, N) 


N) 


which means 


Static Fatigue 


Cov (6,, A) = — 2 * Cov (a, ,6,) 


In = -N In o + In B + (N-2) In S. 


y = a 2 x + b 2 

a 2 = -N , b 2 = In B + (N-2) In S. 
y = In ? f , x = In o a 


given V (a 2 ), V (b 2 ), and Cov (a 2 ,b 2 ) from equations (20) through (22). 
Solving for N in the equation for a 2 gives 


N = -a 2 


Applying the error-propagation rule. 


vvo = (^ y = k <» 2 ) 
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Solving for In B in the equation for b 2 gives 


In B = b 2 - (A-2) In S. = b 2 + (a 2 + 2) In S. 


Applying the error-propagation rule, 

a In B \ 2 


. (£-») „<«,) ♦ (i*±) ribl) + 2 (i^)(4^) cona^) 

+ a In B \ 2 


a In S. 

J 

\2 


-) K(ln§) 


= (In S.) z F (c 2 )+ V (b^ + 2 (In S.) Cov (a 2 ,Z> 2 ) 

+ (a 2 + 2) 2 F (In &) 

Again, assuming the equivalence of F (In 5.) and F (In S), as previously discussed for dynamic fatigue, results 

F (In B) = (In S.) 2 F (a,) + F ((? 2 ) + 2 (In £.) Cov (a 2 ,b 2 ) 

+ (a 2 + 2) 2 i- F (In S.) 

•'i 

Applying the differential expression for covariance, 


m 


Cov (A, In B) = 


d\n B 
dN 

3 In B 


kma - ( d ln B dN , 3 1n B db 2 d j n B 3 ln 

V aiV 3A + a^ 2 dN + ~dN~) F(N) 


3A " In s < 


3A 

aA 


= i 


a in b , a * 2 / </&, . \ 

~3sr = ] - w"(-^ + ins .) (n™' 3 ) 


a In B ,. T a In s. 

JT^. = - ^ 2 ) ' ■' - 


3A 


= 0 


db 


Cov (A, In 5) = -^ F (A) = Cov (3,, A) 


but 


therefore, 


Cov (b 2 ,A) = - Cov (cr 2 ,b ) (Note 4) 


Cov (A, ln B) = - Cov (a 2 ,b 2 ) 
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Note 3: 


db 2 _ db 2 dN 
'dN ~ dN dN 


db 

= s« ■ ln S I 


bb 2 d In B M 2 d In £. 
9 In B dN + d \nS\ dN 

d In S'. 

aS dN~ = 0 


As with dynamic fatigue, using this expression for the total derivative to evaluate the partial derivative results in 
the term -In S. dropping out of equation (A54) in reference 3 as shown in Appendix H. 

Note 4: 

A A 

From applying the error-propagation rule to In B = In B (fl 2 ,b 2 . In S.) and to In B = In B ( N,b 2 , In Sp, 

2 ( TJT K TSf ) C ° V ( “ A) “ S ‘ C ° V ^ 

must equal 

2 C-wf)C-^r) c °" (6 *' N) ’ ( - ln ^ Cm 
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APPENDIX H 


DERIVATION OF DYNAMIC AND STATIC FATIGUE HOMOLOGOUS STRESS 
RATIO ERROR LIMITS NOT ASSUMING A PARTICULAR DISTRIBUTION 


Dynamic Fatigue 


In 


«™ = (^iTr) ta F + JV^T [l " B + , " (W+1,, 

/ 

y - a,x + b , 


a , = 


*3~ ' 3 

1 


3 JV+1 


. K = 


l 


3 N+ 1 


[In fi + In (. N+ 1)] 


’ *-m^ 


y = ln °hd 

given V (a 3 ), V (* 3 ), and Cov (a 3 ,b 3 ) from equations (20) through (22). 
Solving for N in the equation for a 3 gives 


N = — - 1 


Applying the error-propagation rule, 


= (T) V <“3) = (t)^-,) 


Solving for ln B in the expression for b 3 gives 


ln B = ( N+ 1) b r ln (N+ 1) = -^+ ln a 3 


Applying the error-propagation rule, 

,, ,, /91nfix 2 ,, . / d ln fix 2 „ ,, , , - / 9 ln fix / d ln fix _ , , . 

V (ln B) ~ ( da 3 ) F(a 3 )+ ( db, ) V 2 ( da, )( db, ) Cov ( a y b 3> 

b , 


- 3 3 ' '3 

Applying the differential expression for covariance, 


Cov (N, In fi) = 
9 lnfi 


3 ' ' “~3 

- ( - % + + £) + 2 ( - | + i)(i) c<,v m,> 


dN 

9 lnfi 
db , 


d In B 

V(N) = | 

/ 9 ln fi 

dN 

9 ln fi 

db, , 

dN 

l dN 

dN + 

9h 3 

w) ViN > 

( h 

1 \ 

dN 

1 




N+lJ ' 

’ dN 

1 



(A+l) 

db 3 

’ dN = 

db. 

■ dN ~ 

( 1 
\N+ 1 

x*. 

- N \ l) ff 4 ® 1 * 
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db. 


but 


Therefore, 


Note 1: 


Note 2: 


equals 


Cov (N, In B) = V (N) = Cov (b y N) 


Cov (b r N) = -\cov (a 3 ,b 3 ) (Note 2) 


Cov (TV, In B) = - — 2 Cov (a y b 3 ) 


db 3 _ db 3 dN . db 3 a In B db 


dN dN dN + d In B dN 


s d In B \s d \n B \ o, 1 \/ 1 \ 

(-^-)(^ r ) cov w = 2 (-^§ + t)(-z) 


= W + V TV+i K b 3 - N+l ) 


b 3 . 1 \/l 

,2 
3 


2 ( -^-)Hkr) Cov <V0 - 2 W+ »(‘, - i^T ) Cov 

= 2 (y)(6, - a i ) Cov (6 3 ,A'| 


Static Fatigue 


In (iy S^) = - TV In + In 5 


7 = a 4 x + b 4 


a 4 = - TV , T> 4 = In B 


y = In Oy sf) , x = In a HS 


given V (cr 4 ), F (b 4 ), and Cov (a 4 ,b 4 ) from equations (20) through (22). 
Solving for TV in the equation for a 4 gives 


Applying the error-propagation rule, 


Solving for In B in the expression for b 4 gives 


N = -a. 


V (TV) = V(a 4 ) 


In B = b. 


Applying the error-propagation rule, 


F(ln B) = V(b A ) 
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Applying the differential expression for covariance, 


Cov (TV, In B) = d ^j B V (N) = ( 


9 In B 


db. 9TV 


9h, x db. 

) VM-sJvw 


But in this case 


thus, 


9*4 = £^4 
dN dN 

db. 

Cov (TV, In B) = Cov (& 4 ,TV) , 


but 


Cov (b 4 ,TV) = - Cov (a 4 ,b A ) 


thus, 


Cov (TV, In B) = - Cov (a 4 ,& 4 ) 
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APPENDIX I 


DYNAMIC AND STATIC FATIGUE ERROR LIMITS 
ASSUMING WEIBULL DISTRIBUTIONS 

(after reference 3) 


Dynamic Fatigue 

From Appendix F, 

V (A/) =\v(a) = (N + 1) 4 F(a,) 

a \ 

but from equation (A38), page 805, reference 3, assuming Weibull distributions for S and S. and using x f as an in- 
dependent variable results in 


where 


Therefore, 


KJjN-2) 2 

V (0l) " J Q (N+ l) 2 m. R (In 6) 

K d = 1.44 for most glasses and ceramics (Note 1) 

/ = the number of specimens tested at each of J 2 stressing rates <r. for 
a total number of specimens J Q = /,/ 2 

1 Jl . — r 

R (In a) =— S 0 n °r l n °) 
j- l 

1 h 

In a = —r £ In a. 

J 2 7-1 


K d (N-2) 2 (N+ l) 2 
J Q m 2 . R (In a ) 


Note 1: 


K. = m 

a 


2 / N+ 1 \ 

= (w) 


r | 

r i + — < 

f N-2 x-| 

-r 2 | 

[• + 

N-2 I 

L 1 m.y 

<N+ IJi 

m.(N+ 1)J 

r 2 | 

[>■ 

N-2 1 

f m. (N+ 1) J 



+ | In In 2 - m. ^ ) In T [ 1 + m (N+ J } 

T = the gamma function 
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Similarly, from Appendix F, using equations (21) and (22) to substitute for F(&,) and Cov (a,,*,), 
K(lnB) = ( +-L +4lnS,) r(a,) + (l-) [ R (In i) + (tol) 2 ] V (a,) 

' 2 ( _ | + 7 + ^ toS ') (7)i^>'(«,) + (-i--3)V(ln£) 

111 1 W 1 


= "7 b ' + + ln V + ^2 ( ln °) + On a) 2 ] F (a,) 

l a i 

-2 (“j )(— b x + a x + ln S.) InT F (a,) + (— - 3 ) F(ln S .) 

V flj / < 

r i , L , ~ , — , 

~ 2 { _2 (“ + ff i + In S .) 2 + [R (ln o) + (In o) 2 ] 

-|- (- ft, + + !n i) } + (1 _ 3 ) F (ln 5.) 


But from equations (A29) and (A35) in reference 3, 


Expanding the last term in the expression for F (ln B), 


V (ln S) = 

' -V"? 


- 3 K (In S’.) . (!±tL\\ttL\l. n . ,. 2 , * A" « '2 


( A/ 4 - tv ^ 1 

-^ 3 -) (l-3 fll )/?Gn^^-F( fll ) 


Thus, 


F(ln B) = -A {-4 (- ^ + fl, + ln 5. ) 2 + [/? (In a) + ln a 2 ] - i(- b. + a, + In §.) hT* 

(* j ^ Uj d j * ' ( 

2 

+ y 2 ( ) <1 - 3a,) 2 R (In i) } 

VjRlf ,-^+^+1.^ — J r. . , /AT+K .. . % 


Substituting for b 1 and a, gives 


-lna) +/?(lna) [l + 7 2 (-^)(l-3fl 1 ) 2 ]} 


V{a,) 


•{ [- In B - ln (N+ 1) + 3 ln ^ + 1 - In a] + R (In a) [l + J 2 ) ( 1 - -L-)* ] } 
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The factors ( jy _2 ) (l “ N+ 1 ) are approximately 1 , which reduces the expression to equation (A42) in reference 
3 ) 


F(ln B) = 


KAN-2) 2 ( r ^ 1* 

, 2 d n /1 ' < [3 In S, + 1 - In [B (N+ 1)] - In a| 
J Q m 2 R (In a) l L ' J 

+ (1 + J 2 ) R (In a) | 


To evaluate the covariance term, the expression for covariance derived in Appendix F is used, which gives 
Cov (N, In B) = - (N+ 1) (- 5 -) Cov (</,,&,) = ~(N+ l ) 3 Cov (ct^bj 

a i 

Lt 

K^(N+ 2) 2 lna 

Cov ( a.,b .) = -V (a ) In a = - 5 — 5 - — - 

v 1 1 1 J Q (N+ 1 fm 2 R (In a) 


Cov (N, In B) = 


K d (N-l) 2 (N+ 1) In a 
J Q m 2 R (In 6) 


Static Fatigue 

From Appendix F, 


V(N) = V(a.) 


but from equation (A49), reference 3, 


where 


V{a,) = 


K s ( N-2 f 
J 0 m 2 . R (In oj 


= V (N) 


( N -2 ) 2 




[ m. / TV- 2\-i" 

]n In 2 — In r ( ! + — )] 


J 1 = the number of specimens tested at each of J 2 applied stresses, 
<7 ., for a total number of specimens J Q = J l J 2 

R (In a ) = j- £ ( ln % ~ ln a ) 

J 2 j= 1 

1 Jl 

ln = 7 D ln % 

J 2 j—l 
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Similarly, from Appendix F, 


but 


V (In B) = (In §f V (a 2 ) + V (b 2 ) + 2(V In S.) Cov (a 2 ,b 2 ) + (a 2 + 2) 2 F(ln 5.) 
= (In if V(a 2 ) + [R (In a) + hT7 2 ] - 2 In §, hT7 V (a 2 ) 

+ (o 2 + 2) 2 V (In S.) 


K(lni) 


K. 
J.m 2 . 

1 l 


where 


Thus, 


AT. « 1.44 (Note 2) 


K (N-2f 


(N-2) 2 K. 


V (In B) = — / [(In S. - In of + R (In of + - 2 

7 0 /7z/ R (In oj a a m 2 

To evaluate the covariance term, the expression from Appendix F is used, which gives 


Cov (N, In B) = - Cov ( a 2 ,b 2 ) = - (- V (a 7 ) In oj = V (a,) In o n 


K (N-2) 2 In o 
Cov (N, In B) = - s 


J 0 m 2 R (In o a ) 


Note 2: 


7 T (1 + 2/m) - r 2 (l + 1/m) 

K = mj t 1 + [In In 2 - m. In T (1 + 1 /m)] 2 

1 ‘ r 2 (l + 2/m) ' ' 
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